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Abstract 

A few years ago, it had been shown that effects stemming from renormalisation 
group running can be quite large in the scotogenic model, where neutrinos obtain 
their mass only via a 1-loop diagram (or, more generally, in many models in which 
the light neutrino mass is generated via quantum corrections at loop-level). We 
present a new computation of the renornralisation group equations (RGEs) for the 
scotogenic model, thereby updating previous results. We discuss the matching in 
detail, in particular in what regards the different mass spectra possible for the new 
particles involved. We furthermore develop approximate analytical solutions to the 
RGEs for an extensive list of illustrative cases, covering all general tendencies that 
can appear in the model. Comparing them with fully numerical solutions, we give a 
comprehensive discussion of the running in the scotogenic model. Our approach is 
mainly top-down, but we also discuss an attempt to get information on the values 
of the fundamental parameters when inputting the low-energy measured quantities 
in a bottom-up manner. This work serves the basis for a full parameter scan of the 
model, thereby relating its low- and high-energy phenomenology, to fully exploit the 
available information. 
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1 Introduction 


Neutrinos are probably our best handles to detect physics beyond the Standard Model 
(SM) of elementary particle physics, most certainly in the absence of new physics signals 
at the Large Hadron Collider (LHC) and without a clear detection of a Dark Matter (DM) 
particle with well-defined properties. In the last decades, we have experienced tremendous 
successes of experimental neutrino physics. Starting with the first observations of neutrino 
oscillations [TJ, 2|, which proved that at least some light neutrinos must have a non-zero 
mass, we have by now obtained a fairly complete picture of leptonic mixing [3 -T|: all three 


leptonic mixing angles 6^, which describe the mismatch between the neutrino mass and 
flavour bases, have been measured and we have determined two mass square differences 
which constrain the neutrino mass spectrum. Yet, we do not know the neutrino mass 
ordering (i.e., which mass eigenstate is the lightest), we do not know whether neutrinos 
are equal to their own antiparticles (related to an observation of neutrinoless double beta 
decay), and we do not know their absolute mass. The only information we have on the 
latter is that is can at most amount to about 1 eV |8[[9|. 

One possibility to explain the smallness of neutrino masses is by tree-level suppression 
mechanisms, such as the famous type-I seesaw 10-15 . However, a generic problem at 


least with the simplest setting is that it is hardly testable, since the new fields involved 
are generically very heavy. An alternative route to go is to generate light neutrino masses 
at loop-level instead [l6j. In such models, the neutrino mass vanishes exactly at tree-level, 
but it is generated by loop-diagrams as a small correction to the tree-level Lagrangian. 


The easiest such models generate a neutrino mass at 1-loop 17,18 


see 


19 for a general 


classification, although 2-loop 20 -26^ or 3-loop 128-32] mass generation is possible, tooj^] 
The most attractive point of radiative neutrino mass models is that they intrinsically 
connect the phenomenology of light neutrinos to other sectors. The reason is that the 
particles contained in the loop-diagram generating the neutrino mass are typically not 
“invisible”, but rather they have quantum numbers that make them observable. For ex¬ 


ample, the Zee-Babu model 20,21 contains new scalar fields that are electrically charged 
and could thus be detected at colliders 34,35 . Another instance is the Ma-model 


studied in this paper, in which the loop contains particles that could act as DM. In general, 
radiative models prove interesting from a phenomenological point of view because such 
connections to other sectors comprise a handle to distinguish them when combining data 
from different sectors. A generic example is to use input from both low- and high-energy 
experiments and to exploit their complementarity to strongly constrain settings with a 


27 


1 Also for this case a general classification is available, see Ref. 

2 In principle one could go to even more loops, see Ref. 33 for a 4-loop example, however, at that 


point the light neutrino masses would become too small to agree with the experimental bounds. 
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radiative neutrino mass 25 . However, radiative models are also intrinsically interesting 


because of their structure forcing light neutrino masses to depend on many qualitatively 
different model parameters at the same time. This is the key for renormalisation group 
running influencing radiative models in a very non-trivial way. 

Renormalisation group running of such radiative models is still a comparatively young 
held, first mentioned in Ref. 36 . In fact, at least for the so-called scotogenic model 18 .^] 


which is arguably the most simple setting with a 1 -loop neutrino mass and full agreement 
with experiments up to now, the pioneering work on the renormalisation of radiative 
neutrino masses has only been presented in 2012 37 .^] In this reference, the general 
tendencies of the renormalisation of radiative neutrino masses have been worked out: 


• The running can be (very) strong: 

This is easily understood, since loop diagrams depend on products of couplings such 
that, if the two factors ai ;2 in a product a\ x a 2 both receive corrections of the form 
a* —* a* + Aoj, their product receives a correction A(ai x a 2 ) = aq x Aa 2 + a 2 x Acq, 
which can result in enhancements if the terms add up. This fact is trivially reflected 
in the light neutrino mass matrix for the model discussed here, cf. Eq. (Jr]). 

• There is no inconsistency: 

One may ask whether it at all makes sense to work out a 1-loop correction to a 
1-loop diagram leading to a light neutrino mass. However, one has to keep in mind 
that, in fact, what is corrected by the RGEs is not the mass itself but rather the 
Lagrangian of the model. If any diagram of arbitrary order is then computed using 
the RG-improved Lagrangian parameters, one consequently obtains an improved 
result. 


• The origin of lepton number violation becomes clearer: 

A certain set of the model parameters (in the scotogenic model discussed here these 
are the neutrino Yukawa coupling matrix hij , the right-handed (RH) neutrino masses 
Mfc, and the 4-scalar coupling A 5 , to be explained in a moment in Sec. [ 2 ]) drives the 
lepton number violation (LNV) but, if one of them is set to zero, one could in fact 
define a conserved version of the lepton number. Thus, one would expected these 


LNV parameters to be naturally small in the’t Hooft sense 39 , since the symmetry 


of the Lagrangian would be increased if any of these parameters was set to zero. 
This is reflected in the corresponding RGEs, as they only allow for multiplicative 


3 The curious name of this model derives from Ancient Greek and can be roughly translated into 
English as “generated by darkness”, which refers to the DM candidates in this model being part of the 
loop-diagram that generates a non-zero active neutrino mass. 

4 See also the later Ref. 38 for the Zee-Babu model. 
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corrections to the decisive couplings; thus, if one of the LNV parameters is zero 
at any scale, it can never be generated radiatively. This tendency is quite generic 
for radiative neutrino mass models: while running effects can considerably change 
certain observables, they are typically not powerful enough to break lepton number 
in the first place, at least unless the model contains a scalar whose radiatively 
generated vacuum expectation value would do the job. Instead, most radiative 
neutrino mass models contain one or more sectors which intrinsically break lepton 
number, and this breaking is then only translated into the light neutrino sector. 
More technically, this is reflected in the light neutrino mass being proportional to 
all LNV couplings, cf. Eq. Q for the case at hand. 


Furthermore, Ref. 37 has discussed several technical aspects such as how to correctly 
match different effective field theories (EFT) for loop-realisations of the Weinberg opera¬ 
tor 40 or how to correctly integrate out the heavy RH neutrinos present in the scotogenic 


model. 

In this work, we will revisit the renormalisation group running in the scotogenic model 
in a more illustrative way and thereby extend, confirm, and update previous results. This 
manuscript presents a follow-up to 37 which adds several new and important aspects to 
what has been known before. First, while the discussion of the matching in Ref. 37 was 


rather technical, we add a more intuitive picture by dropping some of the more formal 
expressions in favour of graphically matching the different EFTs and directly relating the 
diagrams to the corresponding formula. Second, while 37 relied on a purely numerical 
analysis of the RGEs, we present approximate analytical solutions to the equations for 
several illustrative cases, wherever possible. We relate the approximate analytical results 
to the numerical computations presented, which enables the reader to get an intuitive 
understanding of the behaviour of the quantities involved. And third, by presenting 
results for several limiting cases not discussed in ||37| (all of which are consistent with the 
low-energy neutrino data), we develop an intuitive feeling for the different regimes that 
can appear, depending on the region in parameter space that is considered. 

Note that in this paper we focus on presenting illustrative examples for pedagogi¬ 
cal reasons, which serve to clearly exhibit the various effects running can have in the 
scotogenic model. This means in particular that we compute the renormalisation group 
evolution from a very high to a comparatively low energy scale, as the effects are visible 
more clearly in this case. In reality, however, one may be interested only in a comparison 
between collider and low-energy data sets, in which case the running can be confined to 
this range. Of course, the tools presented also serve for such a case, so that it should 
not be any problem for the reader to restrict our considerations to such a case. Further¬ 
more, different readers may regard some of the examples presented as more “generic” or 
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“natural” than others. While of course such a viewpoint may be very well motivated 
in a concrete setting, we would like to stress that none of the scenarios presented is ex¬ 
perimentally excluded (possibly up to the neutrino oscillation parameters, which is just 
what we want to investigate). This is done purposely, given that we deliberately aim at 
presenting what running in the scotogenic model can do, and what it cannot. A detailed 
phenomenological study where all possible bounds are taken into account in greatest de¬ 
tail is left for future work, in favour of first presenting an illustrative study which allows 
the reader to get a global understanding of the running of radiative neutrino masses. 

This paper is structured as follows: In Sec. [2] we give a brief overview of the scotogenic 
model’s general features and in Sec. [3] we derive analytical equations for the running 
neutrino mass and mixing parameters. In Sec. [4] we discuss the resulting EFTs and 
describe the matching procedure in great detail. We discuss the results we have obtained 
from our analysis in Sec. [5] before concluding in Sec. [6j The 1-loop RGEs can be found 
in App. [A] while Apps. [B] and [C] give some details needed for our computations. 


2 Model Overview 


The scotogenic model is an extension of the SM using several RH Majorana neutrinos 
N\ t 2 ,...; which are SM gauge singlets and can therefore have non-zero masses M 1j2 , ..., 
even before electroweak symmetry breaking (EWSB). At least two of them are needed to 
obtain the two non-zero active neutrino masses which are required to explain the observed 
oscillations of neutrino flavours 41 . We will consider the “next-to-minimal” case with 


three RH neutrinos if not stated otherwise. In addition, there is a second scalar doublet 
r), with SM quantum numbers identical to those of the Higgs doublet. Both types of new 
fields (scalar and fermionic) are odd under a discrete Z 2 symmetry which is imposed in 
order to prevent tree-level neutrino masses and the decay of the lightest Z 2 -odd particle, 
which constitutes a DM candidate if electrically neutral. 


The Lagrangian is given by 18 


l ^ 


c = Csu + 2 Ni$Ni - - (NiMijNf + h.c .) + (D^ (D^rj) - t Lj + h.c.) - R(0, rj ), 

( 1 ) 
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where 77 = ia 2 rj* and a sum over repeated indices is implied. The scalar potential is: 


=m 2 H <ft4> + mlrfr) + y (</>V ) 2 + y (v ] vY + 

+ A 3 (</>V) Wv) + A 4 (0 f r?) 2 + y [(^) (0 + r?) + h.c.] , 


( 2 ) 


where we take A 5 to be real since any phase can be absorbed into 77 . Eqs. ((Tj) and (j2]) 
contain all terms allowed by the symmetries of the model, which in particular do not yield 
a Dirac mass term for the active neutrinos, since 7] cannot develop a vacuum expectation 
value (VEV) due to the Z 2 symmetry. The RH neutrino masses emerge from the matrix 
M upon diagonalising it via held redefinitions. 

Upon EWSB, the Higgs held 0 develops a VEV (0) = (0, v ) T , and one identifies 

the physical scalar helds in 0 = (0,u + yj and 77 = ( 77 + , (? 7 _r + irji) /a/2 ) with the 
following masses: 


m 2 h = 2 AiW 2 = —2m 2 H , 


(3a) 
(3b) 
(3c) 

m 2 = m 2 + v 2 (A 3 + A 4 - A 5 ). (3d) 


m± = m 2 + u 2 A 3 , 

= (^3 + A4 + A5) , 


With these degrees of freedom one obtains for the active neutrino mass matrix from the 
diagram shown in Fig. [l] 




M k h ki h kj j m\ / mji 

i0 S I ji ,r2 


k =1 


32vr 2 


777 


\-Ml V Ml 


m, 


777 2 - ilff 


1 / m I 

° s m? 


E M k h ki h k j 
32n 2 

k=1 


g(M k , 771 r, 777/), 


(4) 


which reduces to the expression: 


(MJij = -J2 y2 J^y hkiM k l f{M k) m 0 )h kj , 


k =1 


(5) 


if the two neutral scalars are nearly degenerate in mass, i.e., if A 5 <C 1. This choice 
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Figure 1: Light neutrino mass generation in the scotogenic model via the exchange of r/° 
and iVi )2) .... 


is motivated by the fact that, if A 5 = 0, we can define a global U( 1) lepton number 
symmetry such that A 5 is naturally small in the’t Hooft sense [39]. Here we have defined 
?71 q = {m\ + rrij)/ 2 , and the function / is given by: 


f(M k ,m 0 ) 


Ml fm |\ Ml 

( M k - mg ) 2 ° g \Ml) + Ml - ml 


( 6 ) 


Note that our expression for J\A. V differs from the original one reported in [l 8 | by a factor 
of 2. This factor was missed in the original reference, but it is required due to rescaling 
the real scalar fields by a factor of 1/\/2 in order for them to be canonically normalised. 
Also, we wish to remark that it is not straightforward to define the functions f or g in 
a basis where the RH neutrino mass matrix M is not diagonal. In order to obtain the 
active neutrino masses and mixing angles for any given energy scale, we need to either 
use numerical tools, or make further assumptions that lead to simplifications of these 
functions. 


3 Analytical Formulae 


In the following derivation of analytical RGEs, we will closely follow previous analyses 
performed e.g. for the seesaw type-I 42,43 , type-II 44,45 , type-III |46|, the inverse 


seesaw 


47 , or for the Weinberg operator 48-50 . The general problem for the case of 


a radiative model is that, in calculating the active neutrino mass matrix, we are relying 
on a perturbative method which in turn requires the RH Majorana mass matrix M to be 
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diagonal (and real). However, owing to the renormalisation group flow, if this is true at 
one scale it need not be true at another. Even if we achieve to define the loop-function 
for matrix valued arguments, the problem becomes manifest when we try to compute the 
derivative of the mass matrix, the latter being schematically written as h T g(M, uir, mj)h: 
a chain rule for a matrix-valued function g is needed, however it cannot be found in 
general unless [. M,M'] = 0, where M' = ^ and t = log We can therefore only 

give analytic expressions in cases where either the loop-function g takes a simple form 
whose derivative can be handled or where its running may be ignored completely]^ 

Let us make one comment on the values of the RH neutrinos before we discuss the 
different limiting cases. In the absence of a concrete mass generation mechanism for 
the RH neutrino masses, their values can in principle be arbitrary. While historically, 
having S'O(IO) Grand Unified Theories (GUTs) in mind 51 , RH neutrino masses have 


been assumed to be much larger than the electroweak scale, there is in fact no physics 
reason for this. For one thing, experiments have not found any sign of a GUT, making 
the argument for huge RH neutrino masses considerably weaker]^] On the contrary, the 
argument based on’t Hooft naturalness, see Sec.[lJ in fact suggest very light RH neutrinos, 
since taking their masses to zero would increase the symmetry of the Lagrangian. In the 
scotogenic model, RH neutrino masses are typically taken to be somewhere at the GeV 

and indeed there are 


to TeV scale for phenomenological reasons (see, e.g., Refs. 53 55] 
arguments for either choice. We thus take on an independent point of view and simply 
discuss several possibilities for the RH neutrino mass. All the examples we present are 
experimentally allowed, but different readers may consider one or the other scenario to 
be better motivated. However, we would like to stress that this paper comprises no 
phenomenological analysis but rather an illustrative and general study of the running in 
the scotogenic model, i.e., our focus lies most on showing examples which demonstrate 
the different effects of the running, so that we do not assume any a priori relation between 
the inert scalar and RH neutrino masses. A stricter phenomenological parameter scan is 
left for future work. 

Aiming to find a suitable classification, some useful limiting cases beyond the assump¬ 
tion A 5 <g; 1 are: 


5 Here, y is the renormalisation scale and go is some reference scale that is needed to make the logarithm 
dimensionless. 

6 As we will see in Sec. [ 5 ] however, ignoring the running can be rather problematic. 

7 In some sense it was not too strong from the very beginning, though, given that some GUTs, such 
as those based on SU (5) (Fl do not even suggest heavy RH neutrinos 


52 
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Limiting case 

C 

p 

i) ml > 

C\ 5 + Ch — C m 2 

i (ft^) +1 (u f u) 

ii) ml ~ 

C\ 5 + c h 

J (ft'ft) +1 (U f U) 

ii) ml ~ Ml 2 

Ca 5 + Ch — \C m 2 

1 (»h) + | (V?U) 

iii) ml <C M, 2 2j 

c Xl , + c h 

| (A*/!) + 1 (pu) 


Tabic 1: Quantities appearing in the analytical RGEs of neutrino masses and mixing 
angles for the three different mass hierarchies. In case ii), ml ~ Mf 2 , we can express 
the mass matrix in terms of m o or M. Which one is used is indicated by bold letters. In 
all expressions we take ml ~ m 2 q . 


i) scalar mass dominates: ml =t M. v — h T ^h, 

ii) similar masses: ml ~ Mf 2 =>- M. v ~ h T M~ 2 h ~ hCm^h, 

iii) fermion masses dominate: ml -C 2 =)> J\A. V ~ — h T M 1 ^1 — log (fyr)) h. 

In the first case we get a suppression due to the large quantity ml = ,n M m i ~ m.f ; , 
while in the second and third case this happens due to large M. In cases i) and ii), we 
have eliminated the loop-function in favour of a much simpler approximation. Case iii), on 
the other hand, still suffers from the aforementioned problem due to the RH mass matrix 
appearing inside a logarithm. However, since the running is logarithmically suppressed, 
one might be led to the conclusion that it could potentially be negligible. 

We are now in a position to compute the RGE for the active neutrino mass matrix, 
which is generally of the form (see App. 01: 

( 4 ■*) 2 M'„ = CM„ + P t M„ + MvP, ( 7 ) 

where C is a flavour-blind function of t and P a matrix with potentially non-trivial flavour 
structure. While C may exclusively influence the running of the mass eigenvalues, only P 
induces running of the mixing angles. Our findings for C and P for the three different cases 
defined above, assuming ml ~ m^, are summarised in Tab. [lj We use the abbreviations 












( 8 ) 


Cp = (47r) 2 rvr for T e {A 5 ,m 2 } (see App. |a| for the corresponding RGEs) and: 


C h = 2T U ~- (gl + 3g%) , 


which is simply twice the flavour diagonal part of ( 47 r) 2 / 2 —1 /i / [see Eq. (A-2b) . 

Following the methods described in the aforementioned references, which are reviewed 
in App. [Bj we have derived analytical expressions for the neutrino mixing parameters. 
To this end we assume dominant neutrino Yukawa couplings, i.e. P = a^h^h + a e Y}Y e ~ 
a.hh)h and h^h = diag(/?, 2 , /?. 2 , h|), where Y e is the charged lepton Yukawa matrix]^ 

The expressions we obtain are: 


(4tt) Wj — 77ii [C + 2 ah [c 12 c 13 /r 1 + S 23 (c 12 Si 3 h 2 T ^ 12 ^ 3 ) A 

+ C 23 ( S 12^2 + C 12 S 13 ^ 3 )] + «/i s i 3 cos 5 sin(20i 2 ) sin (2 6 23 ){h\ - h 2 )] , 
( 47 r)^m 2 = m 2 [C + 2 ah + 523(532833/12 + cf 2 h-3) + 

+ C 23( C 12^2 + s 12 s 13 /? l)] + «/i s 13 cos <5 siii(20 12 ) sin(26> 2 3) (^3 - ^ 2 )] 5 

and 

( 47 r) 2 m[j = m 3 [C + 2 a h [s 2 3 /t 2 + c 2 3 (s 2 3 /t 2 + c 2 3 h 2 )]] 


(9a) 


(9b) 


(9c) 


for the masses, where we have made use of the standard abbreviations Sij = sin 0 l3 and 
Cij = cos Oij with mixing angles 0 l3 . 

For the mixing angles we find: 




(47t) 2 03 2 = -A sin(26 l 32) 

W% 3 = 


m xe 1 ^ 1 + m 2 e 


i<t >2 I 


2 v Am 21 

sin( 202 3 )(h 2 - hf) 


(j 1 ! C 23^2 S 23^ 3 ) + ^(^13), 


A m § 2 


2 1 12 2 + m 3 | 

C32 m 2 e 02 + m 3 


1 + C 


( 10 a) 

+ o{e 13 ), 

( 10 b) 


and 


8 In case ii), TOq ~ 2 , we only find such a simple form for P if we take h to be real or if we replace 

M _1 in favour of 171q . We find that the latter approximation is in better agreement with our numerical 
treatment. 
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X 


( 10 c) 


(47 t) 2 6» , 13 = ^{h\ - h 2 2 ) sin(26>i 2 ) sin(20 23 ) 


m 3 


(1 + C)Am | 2 


X [mi COS(5 - 01) - (1 + C)^-2COs( 5 - 02) - C m 3COs5] + O(0i 3 ), 


Am?, 


where we have used the abbreviation £ = A J 1 • Note that, since there is no mechanism 
that explains leptonic mixing in this model, we have to impose certain values of the mixing 
angles at the input scale by hand. As we will see in Sec. [5j a non-diagonal RH neutrino 


mass matrix M is one possibility to do so, however, Eqs. (10) tell us that a purely radiative 


generation of all mixing angles is not possible. One may also wonder why there is running 
among the mixing angles even though h)h is diagonal. The reason is simply that the 
diagonal entries in this quantity have to be different such that [M. v , P] ^ 0, which then 
induces the running of the mixing angles. 

For the phases we obtain: 


(4vr) 2 h' = 


an(hl + 2 M ( 0 ) + 0 („ 13)i 


26 


( 11 a) 


13 


4tt 2 0 / 1 = a h 


2 777-1777.2 . , , , \ /, 2 2 ; 2 2 j^2\ , 

Cl2 A Sm ^ 2 _ ^ ~ C23 12 ~~ + 


h\ — hi ( 2 777i777 3 2 

+ cos (26 i 23 ) —-— 2 ~ S 12 q -- Sin 0, + c 12 m 2 m 3 sin 0 2 


(lib) 


A 7J7§ 2 


1 + C 


+ 0{d 13 ), 


and 


4tt 2 02 = a h 


2 777 i 7772 . , , ± \ f j 2 2/2 2 / 2 \ , 

S 19 . o Slll(02 01) (h, C 23 h 2 s 23'^3) + 


12 Am 2 , 


+ cos(26 23 )—.—1 s 


Am | 2 '~ 12 


1 + C 


sin 0! + C 2 2 777 2 777 3 sin 02 


+ 0(0 


13 J 


(lie) 


In the expression for <5', we have abbreviated: 


<5 ( 1} = sin( 20 12 ) sin( 20 23 ) 


777 3 


(1 + C)A 777; 


-X 


32 


x [mi sin(<5 — 0 i) — (1 + ()m 2 sin(<5 — 0 2 ) — Cm 3 sin h] 


( 12 a) 


and 
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5<°> = 


mi m 3 


12 


(1 + C )Am § 2 


sin 


(0i - 25) (h\ - s 2 23 h 2 2 - c 23 hfj + 


2 m 2 m 3 . /. 2 2 ; 2 2 ; 2\ , 

S 12 XZ^ sm(0 2 - 25) (/*! - s 23 h 2 - c 23 h 3 ) + 


' Amf 2 


mim 2 

A 777,9 


21 


sin( 0 2 - A) (ft? - 4 ,ft? - s|,fc|) + 


( 12 b) 


+ COs(2$23 


m 3 {hl-h$) 
A m \ 2 


’12 


m\ sin 0 i 

(1 + cT 


+ q 2 m 2 sin 1 


Note that the equations for the mixing angles and (Majorana) phases have been approx¬ 
imated to first order in #i 3 to get more compact expressions, while the ^-dependence of 


the masses is exact. It is evident from Eq. (11a) that for 6*13 —» 0 we run into trouble, 
because 5 is ill-defined at this point. However, as it was pointed out in (501 


5 can be 

continued to the point = 0 by demanding 5' to remain finite. In practice we can avoid 
this subtlety by exploiting this fact and choosing a very small ( 0 13 ~ 10 -3 ) instead of 
vanishing 0 13 whenever this is needed. 


As we can see from Eqs. (10), the running of the mixing angles 0 23 and 0 \ 3 becomes 


large if the difference (h 2 — h 3 ) is large. For 612 to run significantly, we must instead 


ensure 


that the quantity (h\ 


r 2 h 2 

c 23' i 2 


s| 3 /r 3 ) is sizable. Clearly, the enhancement or 


suppression of the running due to the mixing parameters themselves is similar to previous 
analyses done for different models. Especially if the neutrino masses are hierarchical, i.e. 
mi <C m 2 <C m 3 for normal or m 3 <C mi m 2 for inverted ordering, the running is 
strongly suppressed, while for the case of nearly degenerate masses, mi ~ m 2 ~ m 3 , we 
get large enhancements due to the inverse powers of mass square differences. 


4 Effective Operators 


In order to use the convenient MS renormalisation scheme and at the same time extract 
low energy observables, we must take care of the decoupling of heavy degrees of freedom 


by hand, as otherwise perturbation theory breaks down 56,57|. We have done so taking 


into account effective operators of mass dimension up to d = 5 at 1 -loop level. As is well- 
known, in the SM the only relevant d = 5 operator respecting all SM gauge symmetries 
is the Weinberg operator [401, which generalises in a model with N Higgs doublets to: 


(13) 


£c(r = )4" ) CW' W6 4‘e«i0 (,, ‘‘ + ft.<:. 
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Here, c/^ (k = 1,..., N) is one of the N scalar doublets (our convention is = 0 and 
0 ( 2 ) = rj for the N = 2 scotogentic model) and a sum over repeated indices is implied. 
Note that, in the case of the scotogentic model, we have an additional exact Z 2 symmetry 
imposed on our Lagrangian and therefore the operators hi^ and n vanish exactlyj® 

In the present situation, we are interested in the values of the neutrino-related o 
servables at the electroweak scale, i.e. at the Z-boson mass M z — 91GeV{^] Thus, we 
may encounter several mass thresholds between the input scale, which is chosen to be 
equal to the GUT scale M G ut = 10 16 GeV, and M z . Note that we only use the GUT 
scale as one possible example for a high scale, while we do not assume any GUT-inspired 
structure, and in particular no unification of the gauge couplings. Part of the reason for 
this choice is that we simply want to demonstrate what running effects could possibly 
do in the scotogenic model, and this is illustrated much more easily when running over 
several orders of magnitude in energy. On the other hand, a high scale which is at least 
somewhat close to the GUT scale is physically motivated when taking into account that 
the structure of the neutrino flavour sector may be explained by discrete symmetries, 
which are typically imposed at some high scale close to M G ut 58-60 . As we will show 


in the explicit examples presented in Sec. [5j very simple leptonic mixing patterns at the 
high scale could translate into phenomenologically valid regions when running to a lower 
energy scale, which may be of interest for a significant part of the readership. It should 
be clear, though, that in general running effects are suppressed if the running only ex¬ 
tends over a few orders of magnitude in energy scale. Nevertheless, as already argued in 
Sec. [lj the running in models with a radiative neutrino mass can be very strong. Thus, 
the values of certain observables could be modified even if, e.g., comparing their values 
at collider physics energies to low-energy measurements. We will in fact present several 
examples with strong running. Ultimately, we do in this paper provide the general tools 
to analyse the running, which should enable any inclined reader to make use of them 
when e.g. comparing the constraints on the scotogenic model arising from different data 
sets taken at different energies. 

The thresholds we consider are the RH neutrinos TVi, 2 ,... and the inert scalars that 
emerge from the doublet r]. The scheme of resulting EFTs is depicted in Fig. [2] for 
an example mass hierarchy. At each mass threshold M t or m rj , we must remove the 
corresponding field and its couplings from the theory, and match them to the effective 
operators. 

Note that, in the scotogentic model, the active neutrino mass matrix is generated at 


9 It is understood that, by calling re an operator, we are referring to the operator multiplying it. 
10 Below this scale we expect only little running of the mixing parameters, since the Weinberg operator 
is the only source of light neutrino masses and therefore the values at Mz carry over to lower energies to 
a good approximation. 


12 








Figure 2: Schematics of mass thresholds and the resulting EFTs for the case of three RH 
neutrinos and a mass hierarchy Mi < rn v < M 2 < M 3 among the new particles. Below 
the scalar threshold the remaining RH neutrinos no longer couple to the SM fields, such 
that they are effectively decoupled even above their mass thresholds. 


1-loop level, which conversely means that there is no tree-level expression for it. This is 
different e.g. in the various realisations of the seesaw mechanism where the neutrino mass 
is generated via tree-level diagrams or, equivalently, by diagonalising the (active plus 
sterile) neutrino mass matrix 10-15,61-65 . When such a theory is renormalised, one 


can exploit the fact that counter terms of the broken electroweak phase can be obtained 
from the symmetric phase by simple algebraic relations 57 . The same holds for the f3 


functions and all calculations can be performed in the symmetric phase of the theory, 
which is much simpler. 

However, this is in general not true if we are considering a loop-level neutrino mass. 
This becomes obvious if we translate the above reasoning with Lagrangian parameters 
and counter terms into a diagrammatic language: the broken phase degrees of freedom 
are linear combinations of those before symmetry breaking. A tree-level diagram in the 
broken phase can therefore be obtained by a combination of diagrams in the unbroken 
phase. If loops are involved, however, the degrees of freedom propagating in the loops 
cannot be “combined” in such a simple manner - we have to take into account the physical 
degrees of freedom from the very beginning]^] In short, this means that the matching of 


11 Since the masses of the inert scalars are in turn just a linear combination of Lagrangian parameters, 
we can calculate their running masses in the symmetric phase. 
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couplings has to be carried out in the broken phase of the theory. 


4.1 Integrating out the RH neutrinos 

Suppose that we have k RH neutrinos in the full theory. As we run down to lower energy 
scales, we will encounter the mass threshold /j = M k of the heaviest RH neutrino N kl 
and we integrate it out as discussed in detail below. Thereby, we obtain what we will call 
EFT k. Integrating out the next-to-heaviest held yields the EFT labelled (k — 1), and 
so on. In the general EFT n, where (n — 1) RH neutrinos are left, the active neutrino 
mass matrix takes the following diagrammatic form, where we have contributions from 

the Weinberg operator k )( 11) , the Weinberg-like operator ( £ )(22) , and the remaining RH 

neutrinos 0 


(n) 


Vr/i 


</> <b°) 

Tc A 



VLi (™)(H) V Li VL-j (n)(22) V Li V Lj 

K K 


Vr/i 


W, 


, (14) 

"Li 


where the operators W*' 1 have to be matched at all of the previous [k — n + 1 ) thresholds. 
It might seem as though the Erst diagram gives a tree-level neutrino mass, however, the 
matching of the Weinberg operator will reveal that only loop-suppressed diagrams enter 

the corresponding matching condition, such that no loop-diagrams with insertions of k ' ^ 
have to be considered in what follows. In writing it is understood that the diagrams 
with rjji and 77 / in the loop have to be calculated separately and then summed. 

The corresponding analytic expression is: 


M U 2 (n)(ll) 

Muij= “TP 


1 

32vr 2 


l>)( 22 ) 



K ij 

( n , 

w l0g -JF 

2, rnj\ 
— m I l°g -^2 J 

2 



(n) (n) (n) (n) 

hu Mi g(Mi,m R ,m I ) h ej 


£<n 


(15) 


where g is defined in Eq. Q. Note that the mass matrix in the full theory is finite 


12 We assume ordered Majorana masses, i.e. M\ < M 2 < ■ • • < M n . 
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without renormalisation. This has to be true because there is no counter term available 
to cancel a potential divergence of the diagram in Fig. [l]and only renormalisable vertices 
had been used. More technically, in the broken phase this is realised by cancellations of 
the divergent parts in the diagrams for tjr and rji in the loop, stemming from the factor 
i in the expansion y/2 rf = (?/# + irji). At the same time the explicit dependence of 
the mass matrix in the EFTs on /i is signaling that the expression is not independent of 
the renormalisation scheme which we employ, since we have to renormalise the divergent 
(centre diagram) contribution to the mass matrix (a discussion of this can be found in 
Ref. |66|). Of course, once we compare the values of a physical observable between two 
different scales, the prediction will be the same, no matter which scheme is used. 

In practice, we input Yukawa couplings h t j and a RH mass matrix M at the GUT 
scale and run them down to the first mass threshold /i* = Since the running may 

introduce off-diagonal elements in the Majorana mass matrix, we must then diagonalise M 
and cancel the row and column corresponding to the largest eigenvalue in M to obtain the 

(fc) 

effective mass matrix M. Note that, since M must be symmetric, it can be diagonalised 
by a unitary matrix V according to V T MV = diag(M 1; ... ,M k ). By appropriate field 
redefinitions, we find that the Yukawa matrix in this basis reads h = V T h. Removing the 
row that couples the heaviest (in the Majorana mass eigenbasis) RH neutrino N k to the 


(fc) 

lepton doublet yields the new effective Yukawa matrix h . This procedure is continued 
down to the final threshold and it is identical to that typically applied to the type-I seesaw 
mechanism (see e.g. Ref. 67 for further details). 

Note that in the effective theories the mass matrix may consist of up to three indepen¬ 
dent parts [cf. the three diagrams in Eq. (14)], whose running will in general be different. 
Thereby, the running of the mixing angles may be amplified, even in the case where each 
contribution only has flavour-diagonal running. Such threshold effects have their origin in 
divergences that occur in the effective theory, which are however absent in the full theory 
- an effect that is well-known, see e.g. Ref. 42 . Again, Eq. (14) serves as an example: 


while the right-most diagram is finite, the one left to it - which we obtain by removing 
the heaviest RH neutrino - is divergent. 

We now turn to the matching conditions for the effective operators at the RH neutrino 
thresholds. For simplicity, we only consider the matching at the heaviest mass threshold 
M k - the generalisation to others is straightforward. The matching condition for the 
operator at the first threshold is shown diagrammatically in Fig. [ 3 ] Observe that 
the contributions from the fields cancel exactly, since they do not “feel” that a 

heavier particle has been integrated out. Note also that we have to match the operator 
k^ 22 ' > simultaneously since it appears in the last diagram of Fig. [ 3 } The corresponding 
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Vr/i 
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Vr/i 
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VLi V Lj 


(* 0 ( 11 ) 

K 


"Li V L j JVl. fc_t v Li Vl 3 


(k)( 22) 

hi 


VLi 


Figure 3: Diagrammatic matching of the operator kt 11 ) at the mass threshold //* = Mk(n*)- 
Note that no insertions of have to be considered in loops, since this would yield 
effective 2-loop contributions. 

matching condition is shown in Fig. [4| where canceling diagrams are disregarded. 


Since we are working in perturbation theory, we also need to treat the kS u > as per¬ 
turbative quantities which receive corrections from all loop-levels. This suggests that we 
could write: 



(*)(«) ( fc )ho , (mi) , ^ / 2 \ 

K ij — K ij, 0 + e K ij, 1 + o ( e ) > 

(16) 

where e is some small (loop-suppression) factor, e.g. (47 t) - 2 , and l — 1,2. 

By this approach, the leading O (e°) matching conditions can be calculated from Figs. [3] 
and [4] for vanishing external momenta (no sum over k)\ 


i ,v 2 m ip 
° = *y% 0 and 

(17a) 


(fc)(22) 

0 .^ki^kj ! • ^ij, 0 

2 X l ——- = l —— , 

2M k 2 ’ 

(17b) 

and thus: 

K ijt 0 = 0 and 

(18a) 


(k)( 22) ^ hkihkj 

Kij ’° ~ M k 

(18b) 


At O (e 1 ) we obtain, after applying some simplifications and using the tree-level matching 
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Figure 4: Diagrammatic matching of the operator ft* 22 ) at the mass threshold /i* = Mfc (//*). 
Again no loop diagrams with an insertion of ft* 11 ) have to be considered, since they result 
in effective 2-loop contributions. 


conditions (18b) and (|18a|) for the 1-loop diagrams: 
hkihkj 


(fc)(ii) 

K a, i = 


(fc)(22) 

K ij, 1 — ^2 


v 2 M k 

Mi, 


™R m R 


4 

rrij 


l og 122*-22-i 0 g ™L 

m 2 R - M 2 g M 2 k m 2 - M 2 g M 2 


M 2 

3 ( f(M k ,m R ) + log^ 


m 


i? 


Ml 

/(M fc ,m 7 ) + log—| 


ml 


(19a) 

(19b) 


with the auxiliary function / being given explicitly in Eq. (§. 

Putting all the pieces together and generalising to the n th threshold, we find the 
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matching conditions: 


(n)(ll) 
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(n)( 22) 


fi — M n 


(n+l)(22) 
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fl—Mn 


foni hnj 
M n 


2+ ^ U (f(M n ,m R ) + log^ ) - 


(47r) S 
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i? 


M 2 

f(M n , mj) + log^ 


m. 


/i— M n 


(20b) 


In these expressions it is understood that all quantities (especially M and h) are to be 
evaluated in the n th effective theory. Also, since the operator k^ 22 ' 1 only appears in a 


1-loop graph for the mass matrix, see Eq. (14), it is sufficient to consider Eq. (20b) to 
leading order only. 


4.2 Integrating out the inert scalars 

Finally, we must consider the possibility that we have to integrate out the inert scalars, as 
it is shown in Fig. [HJ Since by Eq. (|3]) their mass scale is set by the value of either rn 2 or 
v 2 , we integrate them out at a common scale /i* = m^(^*) if this threshold is encountered 
between the high input scale and Mz- The resulting EFT carries the label “0” in Fig. [2} 

In this case the matching will have several effects: one is that the scalar sector is 
reduced to the simple SM form with only one Higgs doublet. This also modifies the gauge 
coupling RGEs and thereby all other RGEs, which has to be taken into account in the 
numerical analysis. Also, we remove the operator n^ 22 ' > and all neutrino Yukawa couplings 
hij. Therefore, once the inert scalars are integrated out, the RH neutrinos “decouple” from 
the remaining fields and the active neutrino masses are generated exclusively through the 
Weinberg operator k^ 11 ) (cf. Fig. [ 2 J. 

Assuming that (n — 1) RH neutrinos are still present, the matching condition shown 


18 



















</> <0°) 


</) 


Vr/i 


+ 


Vr/i 

S' 

/ 


+ 


v L j W(U) VLi u Lj (n)(11) v Li v Lj N lt , v Li U Lj („)( 22) zy,* 

^ K, K 


Figure 5: Diagrammatic matching of the operator z^ 11 ) at the mass threshold /i* = 


in Fig. [5] gives: 
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fd — M n 


Once this threshold is encountered, no further thresholds need to be taken into account, 
since the RH neutrinos are decoupled. It is interesting to observe that the contribution of 

__ jji 2 

the last diagram in Fig. 5 to Eq. (21) is proportional to the logarithm log fl 2 /7 - which is 


small. In fact, if we matched at a scale rriR ~ mi instead, this contribution would vanish 
and the matching condition would simplify. However, this would require a modification of 
the numerical determination of the threshold, so we refrain from this alternative approach. 


5 Numerical Analysis & Analytical Results 

For our numerical study we have used the full set of RGEs given in App. [A], neglecting all 
SM Yukawa couplings but that of the top and considering the case of three RH neutrinos. 
As we have verihed, this is a good approximation as long as the running of the mixing 
angles is driven by the neutrino Yukawa couplings rather than the charged lepton Yukawa 


couplings. This can also be seen from Eqs. (10) using the replacement a.hh\ 2 3 ^ othh\ 2 3 + 

The purpose of this section is twofold: on the one hand, we illustrate how the running 
of parameters translates into running of physical observables such as neutrino masses 
and mixing angles. We will see that running effects can be large without fine-tuning 
parameters. On the other hand, we wish to compare the results predicted by our analytical 
formulae to those of a more detailed numerical study. Since the analytic equations are 
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valid assuming 6h 3 to be small, we will impose a bimaximal mixing pattern |68j at the 
high input scale (which we take to be equal to the GUT scale for definiteness) such that 
0 13 = 0° and 6 12 = # 23 = 45°. Due to the potentially large running of the mixing angles 
this can yield the experimentally measured values at the electroweak scale. 

In doing so we need to ensure that we meet all requirements for the analytic formulae 
to be applicable. This means especially that h'h must be diagonal. We may achieve this 
by starting off with a diagonal neutrino Yukawa matrix h = diag (hi, h 2 , h 3 ) at the high 
scale, but leaving the RH mass matrix M arbitrary. We can then transform into a basis 
where M is diagonal. The procedure is the following: 

Let us start with the active neutrino mass matrix, which can be diagonalised as: 

U t M u U = D u = diag (mi, m 2 , m 3 ), (22) 


where U is unitary. I11 a basis where the charged lepton Yukawa matrix is diagonal, it is 
just the PMNS matrix. 

With our expression for the mass matrix in the scotogenic model, Eq. ([!]), we get in 
the limit A5 1: 

D v ~ U T h T (- v 2 -^M~ 1 f(M, m„)^ h U = (U T h T ) A (h U ), (23) 

which we can solve for A. I11 general, this matrix will not be diagonal but, since it is 
symmetric, we can diagonalise it with the help of another unitary matrix V: 

diag(A 1; A 2 , A 3 ) =D a = V T AV = V T (U T h T y 1 D v (hU)~ 1 V 

= (V^h) T ~ 1 JJ*D u f/ t (U t /?.) _1 . (24) 


Thus, if we fix m v , the three active masses, mixing angles, and phases at the GUT scale, 


we may solve Eq. (24) for the Ai^ (and thereby In this basis the Yukawa matrix 

is of the form h = UUliag (hi, h 2 , h 3 ) such that h)h is indeed diagonal (as assumed in our 
analytical estimates). 

This procedure is similar but not identical to the Casas-Ibarra parametrisation 


69 


and it allows us to impose that h'h is diagonal, which is needed for the analytic equations 
to be applicable. If this is not required, we may instead fix the RH mass matrix and use 
the Casas-Ibarra result: 

h = y/Ar^Ry/D^U*, (25) 


where now we assume A (or equivalently M) to be diagonal, while R is a (complex) 
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orthogonal but otherwise arbitrary matrix and U is defined as above. 

For the scalar couplings we find a suitable choice to be: 

(Al A 2 , A 3 , A 4i As)| Mgijt = (-0 .01, 0 . 1 , 10-®, 10-®, 1 <T 9 ), 


( 26 ) 


where at large energy scales we need to violate vacuum stability (Ai < 0) in order to 
reproduce the measured Higgs mass and self-coupling, as it is the case in the SM. It is 
a known shortcoming of inert doublet models that it is not simple to ensure stability of 
the scalar potential up to the Planck scale, see Ref. 70 . Since we focus on illustrating 


the consequences of the running in the scotogenic model for the neutrino sector, we have 
chosen the above values in order to avoid the additional difficulty of having to deal with 
an unsuitably chosen potential as well. This could distract the reader from the actual 
messages we would like to bring across. However, we would like to point out that of course 
the running of the scalar potential may introduce all kinds of additional difficulties when 


trying to find a realistic choice of parameters (the parity problem 71 of the scotogenic 
model being just one example). 


5.1 Dominant scalar mass 

In our first example, we show the results of both analytical and numerical computations 
for the case i), i.e. M%. The Yukawa couplings are of 0(1) and we choose a scalar 

mass parameter m v = 350 GeV at Mgut- The RH neutrino masses are fixed according 
to the above procedure and take values < 100 GeV. Besides, we have chosen Majorana 
phases <p\ — // and </> 2 = while the Dirac phase is 5 = 0 at the high input scale. 


Note that such a scenario can be under tight constraints coming from collider physics 72 


and lepton flavour violating (LFV) processes [54,73 , since many such processes ultimately 
yield lower bounds on the RH neutrino and/or inert scalar masses, which can in particu¬ 
lar be dangerous for RH neutrino masses around the Z- pole. However, given the nature 
of LFV diagrams at both low and high energies (where the outer states are comprised 
of SM leptons while the virtual particles involved are the new scalars and/or fermions 
of the scotogenic model), what is actually constrained is only a combination of certain 
masses and couplings. While in a detailed phenomenological study all the collider and/or 
LFV bounds have to be included, the approach we take here for illustrative purposes is 
to assume the neutrino Yukawa couplings small enough that the bounds are no problem. 

Fig. [6] shows the results for this example. First, let us emphasise that indeed strong 
running effects can be achieved without fine-tuning parameters to artificially small or 
large values. This is true for both mixing angles and masses. Note also that below the 
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(a) mixing angles (b) masses 



(c) mass square differences (absolute) 


(d) mass square differences (relative) 


Figure 6: Running mixing angles and active neutrino masses for heavy scalars [case i) 
in section [3] . Solid lines are numerical results and the dotted lines originate from our 
analytical equations. The coloured areas indicate the experimental 3cr ranges for the 
mixing angles and mass square differences in normal mass ordering. 
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scalar threshold, i.e. in the SM effective theory, the mixing angles no longer run 

The analytical and numerical computations agree well as long as #13 is small. However, 
as we go to smaller scales, 0 i3 must grow to reach its measured value, 6 13 ~ 10°, and the 
results of the analytical treatment will inevitably deviate. Nevertheless, they provide a 
useful tool to estimate the values of masses and angles at the low scale given some mixing 
pattern at a high scale, or vice versa. 

From the example shown in Fig. 6(a) we extract the following mixing angles at the 
low scale (ji = M z )\ 


0i2 = 35.78 c 


0 23 = 51.98 c 


0i3 = 7.85 c 


(27) 


The running of the mass square differences is shown in Fig. 6(c) and we can read off: 


A = 7.5 x 10 0 eV , Amg 2 = 1.0 x 10 d eV . 


(28) 


While the angles are well within their experimental 3cr ranges, as is the solar mass square 


difference A'/n 21 74 , the atmospheric Am 32 is too small. I11 fact, for the case of dominant 
scalar mass and diagonal h'h, we have not found a setting where all parameters fit the 
experimental results. This should, however, not be taken too seriously given the strong 
assumptions on the form of the Yukawa coupling matrix h and the RH neutrino mass 
matrix M. Moreover, Fig. 6 (d)| teaches us that the running of the mass square differences 
is not only large in the full theory, but also in the SM as an effective theory below the 
scalar threshold. This is easily understood from Eqs. @> : the factor C contains gauge 
couplings and the top Yukawa coupling, which are both large and also run significantly. 


Finally, we have used the Casas-Ibarra parametrisation, Eq. (25), to show that we 


can indeed reproduce the experimentally determined values of the mixing angles and 
mass square differences when relaxing our requirement that h)h be diagonal. Again, we 
choose rn v = 350 GeV and and RH neutrino masses (Mi,M 2 ,M 3 ) = (50,100,125) GeV 


at M gut . All mixing angles, see Fig. 7(a), and mass square differences, see Fig. 7(c 


are in agreement with the experimental results. Note that the running of 0 2 3 is reversed 
compared to the previous case, which can be explained by the appearance of off-diagonal 
elements in h)h with respect to the analytical RGE, Eq. ( |10a[ ). We see that the running 
of Am 32 is very strong, which is also nicely visible in Fig. |7(d) where the mass square 
differences are plotted relative to their values at the Z mass. The reason for this is that 


m3 runs to very large values, as can be seen from Fig. 7(b) 
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This would still be approximately the case had we not neglected all SM particles but the top quark: 


while Eqs. (10) tell us that the running will then be induced by lepton Yukawa couplings, they are however 


very small and can only cause very little running (cf. [50 
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(a) mixing angles 


(b) masses 



(c) mass square differences (absolute) (d) mass square differences (relative) 


Figure 7: Running mixing angles and masses if the assumption h)h = diag {h\, h|, h§) is 
dropped. All mixing parameters agree with the experiments (lightly coloured boxes mark 
the 3cr ranges). Note the prominent running of the mass square differences, especially of 
A m§ 2 . 


24 






















(a) hlh diagonal (b) Casas-Ibarra 


Figure 8: Running phases in the case of diagonal (left) and arbitrary (right) h)h. Note 
that the phases can be restricted to the interval [0, 2n) in our parametrisation. 

We wish to conclude this subsection by discussing the running of the phases, which 
is shown in Fig. [8} The left panel of this figure shows the running phases in the case 
where we impose h'h to be diagonal. In this case, in particular the running of S is highly 
pronounced, and it is driven from its starting value 5 = 0 to a value of almost 7r, such 
that the running interpolates completely between two regions of CP conservation, while 
passing through a region of maximal CP violation in between. This behaviour is correctly 
captured by the analytical estimates (in dotted lines). The right panel displays the results 
when using the Casas-Ibarra parametrisation, where the running is much weaker. 

5.2 No hierarchy among scalar and RH neutrino masses 

This limit, which corresponds to case ii) in Sec. [3j is quite difficult to realise since, as we 
have learned from Sec.[3j we need 0(1) Yukawa couplings to achieve large running effects 
and therefore at the same time we need a hierarchy among the RH masses to reproduce 
the desired mass square differences. This means that the RH neutrino masses cannot 
be exactly degenerate and the requirements Mi.2,3 ~ m v can hardly be met all at once. 
For our example shown in Fig. [9j we have chosen m v (Mqut) = 700 GeV, which requires 
together with our 0(1) Yukawa couplings and example phases of 0i = if, </>2 = if’ an d 
6 = 0, RH Majorana masses: 

(M 1; M 2 , M 3 ) = (340,407,851) GeV. (29) 
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(a) mixing angles (b) masses 



(c) mass square differences (absolute) (d) mass square differences (relative) 


Figure 9: Running mixing angles and active neutrino masses for scalar masses comparable 
to those of the RH neutrinos [case ii) in Sec. [3] . Solid lines are numerical results and the 
dotted lines originate from our analytical equations. The coloured areas indicate the 
experimental 3cr ranges for the mixing angles and mass square differences in normal mass 
ordering. 
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Clearly, we cannot expect the approximation to hold very accurately and it turns out 
that we find better agreement if we replace M by m v in the formula for the mass matrix, 
as we had suggested in Sec. [3j This means that we use the results of the second row for 
case ii) in Tab. [TJ 

Glancing at Fig. [9] we see that the tendency and the direction of the running are 
correctly captured by the analytical approximations, however, the running is largely un¬ 
derestimated. This can be explained by two effects. First of all, as mentioned above, 
approximating the mass matrix as we have done cannot be very accurate. And secondly, 
since 0 i 3 is growing faster than before as we lower the renormalisation scale, the assump¬ 
tion of a small #13 breaks down more rapidly. 


This time at /j = Mz we numerically find the mixing angles, see Fig. 9 (a 


012 = 31.36°, 023 = 49.29°, 0 , 3 = 8.85°, 

and mass square differences, see Fig. 9(c)[ 

Am 21 = 7.5 x 1CT 5 eV 2 , Am 32 = 4.6 x 1CT 4 eV 2 . 


(30) 


(31) 


As before, only the solar Am 21 agrees with the value reported by experiments. The 
good agreement between numerics and analytical results for the mass square differences 
is rather surprising, since neither the angles nor the masses themselves suggest such good 
agreement. Most likely this is simply a lucky choice of input parameters. 

Again, we have also found a less restrictive example for similar parameters, using 


the Casas-Ibarra parametrisation, see Fig. 10 which reproduces all experimentally de¬ 
termined oscillation parameters. Here, we input m. q = 240 GeV and (Mi,M 2 ,M 3 ) = 
(250,320,960) GeV at M GUT , as well as phases </> 1 ~ |, (p 2 — ^y, and 5 ~ ^y. Two fea¬ 
tures are worth being highlighted at this point, one of which is that for the first time we 
have a RH neutrino threshold above the scalar threshold, which leads to observable but 
not overly strong threshold effects. Furthermore, two of the light neutrino masses almost 


meet at a scale /1 ~ 10 12 GeV in Fig. 10(b), which drives the running of $12 to extremely 


small values, cf. Eq. (10a). Eventually, the difference between the masses grows again as 


does 0 i 2 - This is a nice example of how degeneracies (here generated by radiative effects, 
rather than by a choice of the input parameters) may strongly drive the running. 


5.3 Dominant RH neutrino masses 

Finally, we discuss the case of large RH masses, however, we do not even attempt to 
compare this scenario to an analytical estimate since it is doomed to fail. This can be 
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(a) mixing angles 


(b) masses 



(c) mass square differences (absolute) (d) mass square differences (relative) 


Figure 10: Running masses and mixing angles in case of no mass hierarchy among the 
new particles if the assumptions on h)h are relaxed. 
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(a) mixing angles 


(b) masses 



(c) mass square differences (absolute) 


(d) mass square differences (relative) 


Figure 11: Running of the mixing parameters for dominant RH masses. Note 
that the scalar mass is “attracted” to the (TeV scale) RH masses via the last 
Eq. (|A-5bl). 


however 
term in 
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easily understood by recalling that we need to assume all quantities appearing inside the 
logarithm not to run. However, even though logarithmically suppressed, the running of 


rn 2 cannot be ignored, since the last term in Eq. (A-5b) makes rn 2 run over many orders of 
magnitude if Mi i2i3 m n and h t j ~ 0(1). Therefore, choosing a small m n at Mgut will 
result in an m v comparable to the Mi . 2i 3 at the electroweak scale - and the approximation 
breaks down. A way out would be to choose m 2 < 0 at the high scale, but this may result 
in breaking the Z 2 symmetry at high scales and hence the expression for the mass matrix 


being meaningless 71 


Fig. 11 shows an example for such a scenario. Here, we have chosen m„ = 240 GeV, 


Mi ~ 500 GeV, and M 2j3 ~ 1 TeV. The phases are input as <pi ~ S ~ ^ and 0 2 ~ 
The dashed grey lines indicate the thresholds as labelled in the plot. As expected, the 
scalar mass m r/ is very sensitive to the TeV scale RH masses, and at yU* = m^/x*) it is of 
the same order (it even exceeds the smallest mass, which is why the threshold Mi is not 


plotted). Note the strong running of the mass square differences in Figs. 11(c) and 11(d) 


In spite of the aforementioned obstacles, the example shown reproduces the oscillation 
parameters correctly. It is also worth noting that in this case the threshold effects are 
quite significant, as discussed in Sec. [4| 


5.4 Inverted mass ordering 

The previous analyses are of course equally applicable to the case of inverted mass or¬ 
dering. It is, however, more difficult in our approach to find results in agreement with 
experiment because, since m 3 is the smallest mass, the running of #13 is damped [cf. 
Eq. (10a)]. Therefore, generating a non-zero #13 exclusively via radiative corrections, as 


we have chosen to do so far, is very difficult. We show an example of this in Fig. 12 
where we have a similar setting as in the case of dominant RH neutrino masses. We 
input m v = 120 GeV and (Mi, M 2 , M 3 ) = (450, 850, 900) GeV at Mgut, and as before 
we fold that the scalar mass is attracted to the scale of the RH neutrino masses at lower 
energy scales. All phases are zero at the input scale in order to achieve the largest possible 
value of 6 13 (see the next subsection for a discussion of the effects the phases have on the 
running). 


The qualitative features of Fig. 12 are easily understood. We start with a rather small 
Am|i rs./ 10 4 eV 2 , which grows very fast and then slowly decreases to its value at M z 
[Fig |12(cj] . This fast growth is driving the extreme running of 0 12 in Fig. 12(a) at scales 
just below Mgut- Further lowering the scale /x has virtually no effect on 0 12 , because 
A mh is larger than before and the neutrino mass scale is monotonously decreasing, thus 
suppressing the running of all the mixing angles. 
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(a) mixing angles 



(c) mass square differences (absolute) 



(b) masses 



(d) mass square differences (relative) 


Figure 12: Running of the mixing parameters for inverted mass ordering. Since m 3 
is the lightest among the masses, the running of #13 is not strong enough to reach the 
experimentally preferred region. 
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(a) 


» 1,2 


= 0 


(b) ^ = 6 = 0 


Figure 13: Running mixing angles in a bottom-up approach in normal ordering, where 
the values at the Z-boson mass have been hxed to their best-fit values 
(j>i,(j> 2 , and 5 are hxed at this scale, too, as indicated in the plots. 


74 . The phases 


5.5 An alternative approach 

When studying the scotogenic model’s parameter space, it may be more convenient to use 
a bottom-up approach as opposed to the so far employed top-down approach. The reason 
for this is simply that in principle all model parameters should be considered as input at 
the high scale in the latter case. 14 Obviously, this increases the size and dimensionality of 


the parameter space, which can be very costly in a phenomenological study. In contrast, a 
bottom-up approach allows one to fix all low-energy observables to their measured values. 
To ease the computation and remove some of the ambiguities of such an approach, we 
only consider the case where the three RH neutrinos and the inert scalars are almost of 


the same mass, slightly below 1 TeV (vertical dashed line in Fig. 13). A convenient choice 
for the scalar couplings is: 


(A 2 , As, A„, As) |„ =Mi , = (0.1, 0.1, 0.2, IQ- 8 ), 


(32) 


which guarantees perturbativity, an intact Z 2 symmetry, and the (meta-)stability of the 
scalar potential up to at least 10 16 GeV. 


I 11 Fig. 13 the results of such an approach are shown. Since we expect no qualitatively 


14 In practice however it turns out that it is sufficient to treat the Higgs mass as a free input parameter 
at the high scale and run the remaining couplings from the low to the high scale and back. The changes 
due to this approach lie within the experimental uncertainties in almost all cases. 
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new behaviour, we would like to use this to illustrate the role of the CP phases 0i,. 


and S. From Eqs. (10) we learn that a non-zero S can only influence 6 \ 3 directly, which is 


confirmed by Fig. 13(a) while 0 2 3 is largely independent of the choice of S, 0 12 is modified 


only indirectly via changes in $ 13 , which enter Eq. (10a) through higher order corrections. 


Non-zero Majorana phases 0 i t2 can cause a suppression of the running, as we can see 


from Eqs. (10a) and (10b). For example, the running of 612 is suppressed compared to 


the case of vanishing or equal Majorana phases if we have |0i — 02 1 = tt, because in that 
case the decisive factor becomes: 


rriie 1 ^ 1 + m 2 e l< ^ 2 1 
A m 2 2l 


—>■ 


\mi — m 2 \ 2 \rrii + m 2 


A?t4 


< 


AmJj 


(33) 


As an example, we direct the reader’s attention to Fig. 13(b)| which illustrates this for 
the case where only 02 can take non-zero values at the low scale /i = Mz■ As we vary 
02 from 0 to 7 T, we achieve a less pronounced running for all mixing angles and thereby 
(accidentally) achieve a unification of the mixing angles #12 and $23 at the GUT scale, 
for 02 = |. For 0 2 = n we see that the running is essentially turned off due to the 
aforementioned suppression mechanism. 


6 Conclusions 

We have studied the scotogenic model, a particularly simple setting that generates active 
neutrino masses at 1-loop level and exhibits several Dark Matter candidates. In doing 
so, we have updated several previous results, namely the formula for the light neutrino 
mass matrix, the matching conditions at the particle thresholds, and the 1-loop RGEs 
of the model. We have derived analytical equations that allowed us to identify regions 
in parameter space which we suspected to show significant running behaviour. This was 
verified by comparing the analytical estimates to a numerical treatment. 

The point of this paper is not so much to perform a very detailed phenomenological 
study, but we rather aimed at an illustration of the effects the renormalisation group 
running can have on the scotogenic model, and how to understand them using approxi¬ 
mate solutions of the evolution equations. Following this approach, we have intrinsically 
disregarded some points which could be very crucial in a realistic analysis. These include 
motivating the mass patterns used, explaining the leptonic mixing pattern by a concrete 
symmetry, or presenting an explicit justihcation for the scales used in the various exam¬ 
ples. We nevertheless consider our approach to be valuable, because it will enable us 
- and hopefully any reader - to make use of our general results in a concrete study. 

The key results obtained in this work, apart from clearly illustrating how strong run- 
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ning effects can be in the scotogenic model, are the updated RGEs, their analytical 
approximate solutions (none of which have been known previously), and in particular the 
explicit discussion of all qualitatively different scenarios that can appear. Our work lays 
the foundation for a detailed study aiming at a determination of the full allowed param¬ 
eter space in the scotogenic model. Given that such an endeavour can only be done in a 
purely numerical manner, having a picture of all running effects that could possibly ap¬ 
pear as well as some limiting cases, which can be viewed as “cornerstones” for a numerical 
computation, will be highly useful. 

The scotogenic model is the prime example for a setting with a radiative neutrino 
mass and it has been studied from all sides, from low-energy measurements over neutrino 
physics and Dark Matter constraints to collider phenomenology. However, there exists 
up to now no comprehensive work, which really tries to put together all constraints and 
work out which regions in the parameter space survive. A few attempts exists which e.g. 
confront collider data with Dark Matter signatures, but renormalisation group running is 
completely disregarded in such cases. This study will hopefully contribute to close all these 
gaps by our equations providing the basic tools for a fully comprehensive phenomenological 
analysis of the scotogenic model which can take into account the effects of having data 
available at several different energy scales. 
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A Renormalisation group equations 


The 1-loop RGEs for the scotogenic model have first been computed in Ref. 37 . We 


have re-derived those equations needed for the purpose of this paper, and have in passing 
taken the opportunity to update part of the earlier results. 

For convenience, we define the differential operator T> = The 1-loop RGEs 

for the gauge couplings are those of a generic two Higgs doublet model (THDM) [75]: 


Vgi = hgl 


no sunn 


(A-l) 


with b = (7, —3, —7). 

The quark sector of the scotogenic model is the same as that of the SM, such that the 
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corresponding RGEs do not change P’1 The RGEs for the leptonic Yukawa couplings are: 


VY e = Y. | IyJy c + ifctfc + T - jgl - | 9 || , (A-2a) 

Vh = h | + \v}Y e + T„ - , (A-2b) 


where T u = Tr (h)ti) and T = Tr ^1^17 + 31^ + 3Y^Y^j. For the Majorana mass 

matrix, one finds |37,67|: 

VM= {{hh ] ) M + M (kb))*}. (A-3) 


For the quartic scalar couplings, we find the RGEs for a Z 2 symmetric THDM 


76 : 


V\i — 12A| + 4Ag + 4A 3 A 4 + 2A4 + 2Ag + - (g 4 + 2glg% + 3g 2) 

~~ 3Ai [g\ + 3^2) + 4AiT — 4T 4 , 

V\ 2 = I2A2 + 4A3 + 4A 3 A 4 + 2A4 + 2Ag + - (<? 4 + 2glgl + 3^0) 

— 3A 2 (g 2 + 3^2) + 4A 2 T u — 4T 4l/ , 

VX 3 — 2 (Ai + A 2 ) (3A 3 + A 4 ) + 4A3 + 2A4 + 2A5 + - (g± — 2g\g.\ + 3^2) 

~~ 3A 3 {(j\ + 3^2) + 2A 3 (T + T v ) — 4 T ve , 

7AA 4 = 2 (Ai + A 2 ) A 4 + 8A 3 A 4 + 4A 4 + 8A 3 + 3g±g 2 

— 3A 4 [g 2 + 3^2) + 2A 4 (T + T y ) + 4 T ue , 

T>\ 3 = A 5 [2 (Ai + A 2 ) + 8A 3 + 12A 4 — 3 (57 + 3^2) + 2 (T + T„)], 


(A-4a) 


(A-4b) 


(A-4c) 


(A-4d) 


(A-4e) 


where we have used the abbreviations T 4 = Tr ^Yj*Y^,Y4Y e + 3Y£Y u Y£Y u + 3Y f jY d Y^ Ydj , 
T 4i , = Tr (ftt/i/it/i) and T„ e = Tr (tfhYjY e ). 


15 Note the implicit changes in g-\ .2, though, by virtue of Eq. (A-l). 
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The scalar mass parameters obey the following RGEs: 


Vm 2 H = 6Ai??% + 2 (2A 3 + A 4 ) rn 2 + m 2 H 


2 T-- ( 9 ? + 392 ) 


(A-5a) 


Vm 2 = 6A 2 m 2 + 2 (2A 3 + A 4 ) m 2 H + m 2 r 


2 T„ ~■ | ( 9 ? + 3 9 I) 


- 4 £ 

(A-5b) 


2=1 


where the last term in Eq. (A-5b) is nothing but a trace and thereby invariant under the 


transformation that diagonalises M, such that we do not have to perform this diagonali- 
sation explicitly. 

Finally, we obtain for the effective operators: 


(n)(ll) 1 ( W'H 


(n)(ll) 1 (n)(ll) ( W'W 


M h h - 3YX ’ ( h h - 3Y e f Ye 1 + 


2 \ ( n )( n ) . 0 , ( n )( 22 ) 


+ (2 T + 2Ai - 3 gl) k + 2A 5 


K 


(A-6) 


vT 2) = \ (A ft + n'u) ( k ,( 22) + i ( K ,(22, f ft' ^ + n'u 1 + 


(n)(22) (n)(ll) 


+ (2 T„ + 2A 2 - 3 gl) k + 2A 5 


K 


(A-7) 


B Derivation of the Analytical Formulae 

In this appendix, we review the derivation of analytical equations for the running mixing 
angles and masses. As mentioned in Sec. [3j this is analogous to previous results found 
e.g. in 42-50 . 

To arrive at an analytical expression, we define t = log (in the following a prime 
denotes the derivative w.r.t. t). The light neutrino mass matrix RGE is most generally 
given by: 

(470^ = P t M„ + M U P + CM V , (B-l) 

with C a factor that has no flavour structure, and P a Hermitian matrix. 

Since we generate Majorana neutrinos, the light neutrino mass matrix must be sym- 
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metric and can therefore be diagonalised by a unitary matrix U according to: 


U T M U U = D u = diag(mi,m 2 ,m 3 ), 


(B-2) 


where we may right-multiply U with a diagonal matrix of phases to render the m; real 
and positive. 

Thus, we have: 

D' v = U' t M u U + U t M' u U + U t M u U'. (B-3) 

We can always dehne an anti-Hermitian matrix T such that U' = UT , and the remaining 
task is to determine the components of T in terms of known quantities from the equation: 


D' = D V T -T*D u T 


(4tt) 


CD U + P t D u + D V P 


(B-4) 


In this last expression, we have dehned P = U^PU. 

Since we can choose the Majorana phases in U such that m; e M + for all t, we obtain 
the mass RGEs from the real parts of the diagonal entries in (B-4) (no sum over i): 


rrii 


Til: = 


(4vr) 


Re(C) + 2Re(Pjj) 


and since the imaginary part of the diagonals must vanish: 

1 


2Im(Tjj) = 


(4tt) 


Im(C) + 2Im(P«) 


(B-5) 


(B-6) 


The real part of T it vanishes due to the anti-Hermiticity condition. 

The off-diagonal elements of T can be constructed from (again, no sum over i or j): 


0 = D' vii = 171^ - T-.m- + 


UIJ I (A^\2 


niiPij + nijPu 


3* 3* 


(B-7) 


Taking the real an imaginary parts of this equation, we obtain for a Hermitian P: 

1 m-i + rrij 


Re(T^) — 
Im(Ty) = 


(47t) 2 rrii — rrij 

1 rrii ~ m j 


Re(Pij), 
Im (Pij). 


(B-8a) 


(B-8b) 


(47t) 2 rrii + rrij 

The equations resulting from U' = UT can now be used to extract the running of the 
mixing angles and phases, while the running masses are obtained from Eq. (B-5). 
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C Extraction of Mixing Angles & Phases 


The PMNS matrix can be parametrised in the following standard way 



U 


/ e iSe 

0 

0 ^ 

1 J 

/ g-* 0 l / 2 

0 

°\ 

0 


0 

M 

0 

g—*02/2 

0 

V 0 

0 

e l5r j 


V 0 

0 

1 


with 


/ C12C13 S12C13 si 3 e~ lS \ 

v = -C23S12 - S23Sl3Cl2e J<5 C23C12 - S 23 Sl3Sl2e l5 S23C13 

V S23S12 - C 2 3Si3Ci2e*' 5 -S 23 Ci2 - C 23 S 13 Si 2 e lS C23C13 / 


(C-l) 


(C-2) 


where c t j = cos % and = sin Oij. 

I 11 our numerical treatment, we diagonalise the active neutrino mass matrix according 
to Eq. (B-2) for which we have adapted the algorithm described in 77 to arrive at 


positive and real irii y2)3 . We can then extract the mixing angles and phases according to 


the following relations 42 78 


@12 — 


arctan 


U 12 

Uu 


if U n ± 0, 
else, 


6*23 = 


arctan 


U 23 

U 33 


if U 33 ± 0, 
else, 


6*13 = arcsin \U V3 \, 
— arg (U 23 ), 

S T = arg (C/33), 


6 


arg 


• U^UijUjiUh 

C12CJ3C23S13 


+ C12C23S13 ^ 


S12S23 


for i ^ j, 


(C-3a) 


(C-3b) 


(C-3c) 

(C-3d) 

(C-3e) 


(C-3f) 


S e =arg (e lS U 13 ) , (C-3g) 

0 ! —2 arg (e^t/n) , (C-3h) 
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fa =2 arg (e iS «U* 12 ). 


(C-3i) 


Note that the phases S e , 5^, and S T are unphysical in the sense that they can be absorbed 
into the three left-handed lepton doublets 0 li , 2 ,3 and do not appear in the analytical 
expressions for the running mixing angles [Eqs. (10)] or masses [Eqs. 
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